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Abstract 

We have investigated QCD with two flavors of degenerate fermions using a Symanzik- 
improved lattice action for both the gauge and fermion actions. Our study focuses on the 
deconfinement transition on an Nt — 4 lattice. Having located the thermal transition, we per- 
formed zero temperature simulations nearby in order to compute hadronic masses and the static 
quark potential. We find that the present action reduces lattice artifacts present in thermody- 
namics with the standard Wilson (gauge and fermion) actions. However, it does not bring studies 
with Wilson-type quarks to the same level as those using the Kogut-Susskind formulation. 

C0L0-HEP-381, FSU-SCRI-97-17 
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1 Introduction 



An understanding of the high temperature behavior of QCD is desirable in addressing problems 
such as heavy ion collisions and the evolution of the early universe. It is believed that, at a tem- 
perature between 140 — 200 MeV (where pions are produced copiously) , hadronic matter undergoes 
a transition to a plasma of quarks and gluons. This phenomenon is intrinsically nonperturbative, 
and Monte Carlo lattice simulation provides the best theoretical tool with which to study it. 

Most lattice studies have used Kogut-Susskind (KS) fermions because of their exact U(l) chiral 
symmetry at finite lattice spacing. The full SU(2) chiral symmetry is recovered in the continuum 
limit. Although in contrast Wilson fermions explicitly break chiral symmetry, the continuum limit 
of the different discretizations is expected to be the same. One advantage of simulating with two 
flavors of Wilson quarks versus two flavors of KS fermions is that the updating algorithm is exact 
in the former case but has finite time step errors in the latter. The price to be paid is that the 
explicit chiral symmetry breaking gives rise to an additive mass renormalization; thus, the location 
of the chiral limit is not known a priori. 

Even more troublesome for dynamical Wilson fermions is the presence of lattice artifacts which 
qualitatively affect physics at large lattice spacing. In Refs. p| and [Q] it was found that the 
deconfinement transition becomes very steep for intermediate values of the hopping parameter, k. 
In fact, on an Nt = 6 lattice the transition appears to be first order for a range of intermediate 
hopping parameters and smooth otherwise. 

The lattice community has worked very hard recently to construct actions which have fewer 
lattice artifacts than the standard discretizations of the continuum action. One philosophy is to add 
operators to the action which cancel 0{a n ) terms in the Taylor expansion of spectral observables. 
It is plausible that an action which converges to the continuum action faster in the a — > limit 
would be free of the artificial first order behavior. We adopt this improvement program, attributed 
to Symanzik, in the present work. 

An alternative approach is to search for an action which lies on or near the renormalized 
trajectory of some renormalization group transformation. Since all irrelevant couplings are zero 
along the renormalized trajectory, actions there have no scaling violations, i.e. are quantum perfect. 
Such an action was approximated by Iwasaki § and has being used to study QCD thermodynamics 
with two flavors of unimproved Wilson fermions [Q] . Although it is still an open question to what 
extent this action lies on a renormalized trajectory, the results of Ref. H show improvement over 
standard Wilson thermodynamics. 

In this paper we report on our simulation of finite temperature lattice QCD with two flavors of 
O(a) Symanzik-improved fermions and 0(a 2 ) Symanzik-improved glue. We describe the action we 
used in Section ^. In Section || we give the details of our simulations, and we present our results in 
Section ^j. Finally, we give our conclusions in Section [|. 

2 Action 

When one expands a lattice operator in a Taylor series about zero lattice spacing a, one recovers 
its relevant (or marginal) continuum operator plus higher dimensional irrelevant operators propor- 
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Figure 1: The three Wilson loops in the one- loop Symanzik improved gauge action that we used. 



tional to powers of a. Symanzik suggested that by selecting a favorable combination of lattice 
operators in the lattice action, one might have cancellations of the irrelevant operators up to some 
order in the lattice spacing Q. Liischer and Weisz have applied this philosophy to S\J(N) gauge 
theories. They imposed an on-shell improvement condition whereby discretization errors are elim- 
inated order- by-order in a from physical observables and constructed an 0(a 2 ) improved gauge 
action ||. Furthermore, they computed the coefficients of the operators in the action through 
one-loop order in lattice perturbation theory {?]]. This improvement condition does not provide a 
unique action. The choice which is the most efficient in terms of computational effort adds a 1 x 2 
rectangle and a 6-link "twisted" loop to the Wilson plaquette action. (See Fig. [I].) 

It is well-known now that lattice perturbation theory in the bare coupling go is not trustworthy. 
The a 2 in the vertex of the tadpole graph is cancelled by the ultraviolet divergence of the gluon 
loop. Therefore, hidden in the higher order terms of the expansion in a are tadpole graphs which 
give an effective o° Y^ n c n9 2n contribution. A standard way to deal with this problem is to define 
a mean link no and replace — > U^/uq [||, [|. This introduces a "boosted" coupling constant 

9o/ u i 

Here we combine these two ideas, Symanzik improvement of the action and "tadpole improve- 
ment" of lattice perturbation theory. Our gauge action for this work is as derived in Ref. 
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The subscripts "plaq" , "rect" , and "twist" refer to the lxl plaquette, the planar 1x2 rectangle, 
and the "x, y, z, —x, —y, —z" loop, respectively. Following Ref. || we have chosen to define the 
mean link uq through 

u = ReTr (U plaq )) X '\ (4) 
and the strong coupling constant is defined through the perturbative expansion of the plaquette pi 



2 In (i Re Tr (C/ plaq ) 



T_ 

4tt ~ 3.06839 



(5) 



In the Monte Carlo simulations, we tune uq in the action to be consistent with the fourth root of 
the average plaquette. This procedure is discussed in more depth in Section ||[ 



Figure 2: The clover term in the Sheikholeslami-Wohlert action. 



The Wilson fermion action has errors of 0(a). The Symanzik improvement program can be 
extended to improving this action too. In the a — > limit, the 0(a) term in the action is propor- 
tional to ipD 2 tp and can be removed by adding a next-nearest neighbor operator to the action [12], 
or, after an isospectral transformation of the fermion fields, by adding a magnetic interaction [13]. 
Thus, the tree-level Symanzik improved action is 



S 



f 



"0 x n<v 



ijj(x) ia lxu F IJyU ip(x) 



where Syy is the usual Wilson fermion action, and 



(6) 



(7) 



f^ u is the clover-shaped combination of links (see Figure |2|). As with the gauge action, the links 
here are also tadpole improved. Note that one factor of uq is absorbed into the hopping parameter. 
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Both these gauge and fermion actions have been used widely, e.g. in studies of finite temperature 
SU(3) HH, 15] and quenched spectroscopy [16, |l7|, |j~8|| . At least one group is in the progress of using 
this action to calculate spectroscopy on unquenched configurations [p~9|] . Therefore, we believe our 
choice of action to be well-justified and useful for comparison to other work. 

Of course, further progress is being made in refining the Symanzik improvement program. One 
can attempt to set the coefficients of the higher dimension operators nonperturbatively by de- 
manding that, for example, Ward identities be satisfied up to some order in the lattice spacing pp| . 
Also, fermion actions which are constructed to have errors of C(a 3 ) to C(a 4 ) are currently being 
tested HI |2|. 



3 Simulation details 

Our finite temperature simulations were on an 8 3 x 4 lattice. At fixed (= 6.4, 6.6, 6.8, 7.0, 7.2, 
7.3, and 7.4) we varied k in small increments across the crossover. Our microcanonical time step 
was such that the acceptance rate was between 60-80%; typically At = 0.03, but for the stronger 
gauge coupling (lighter quark mass) runs At = 0.01. We accumulated over 1000 trajectories at 
points close to the crossover. For the updating, we used the standard hybrid molecular dynam- 
ics (HMD) algorithm [23] followed by a Monte Carlo accept/reject step (hybrid Monte Carlo, or 



HMC) []2j]. The calculation of the "clover" contribution to the HMD equations of motion is te- 
dious but straightforward.^] The even-odd preconditioning technique developed for standard Wilson 



fermions [26] is also implemented for the improved Wilson fermion action (27 

Our companion zero temperature runs were performed on an 8 3 x 16 lattice at five (/3, k) points 
along the crossover, as well as at five other points neighboring the crossover line. Some of these 
runs were extended in order to be able to extract the heavy quark potential from Wilson loop 
expectation values. As with the finite temperature runs, we tried to maintain a Monte Carlo 
acceptance rate around 60-80%, so our time step varied from At = 0.005 to 0.03. Measurements of 
hadron correlators and Wilson loops were taken every 10 trajectories. 

The large majority of our simulations were performed on the IBM SP2 at the Cornell Theory 
Center, and two finite temperature simulations were run on a cluster of IBM RS/6000's at the 
Supercomputer Computations Research Institute of Florida State University. 

We used the conjugate gradient (CG) matrix inversion algorithm to compute (M^M) -1 with a 
maximum residue of 10 -6 during the HMC updating. For the spectroscopy calculations, where we 
wish to invert M, we found the stablized biconjugate gradient (biCGstab) algorithm to be twice as 



efficient as CG [28 



We tuned uq so that it agreed with the fourth root of the space-like plaquettes that we measured. 
It might have been preferable to do this on T = configurations, however, the heavy cost of 
repeatedly equilibrating a Nt = 16 lattice forced us to perform this tuning procedure on the Nf = 4 
configurations. In the next section we will show that the difference in the two ways of tuning uq 
is small. This tuning procedure would be dangerous if the system underwent a first order phase 
transition, but, we will also show that the plaquette varies smoothly across the transition. Let 

independent of our work, calculations of the HMD equations of motion for this improved Wilson fermion action 
have been published in |]25|. 
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us remark that this tuning procedure is a prescription, uo may be defined in a number of ways 
since it is an estimate of the higher order tadpole contributions to perturbation theory calculations. 
Therefore, while one might argue that tuning uq strictly on a zero temperature lattice would better 
estimate the tadpole contributions, our method is well-defined and self-consistent. 

4 Results 

4.1 Thermodynamics 

The first task was to locate the thermal crossover line kt{(3). To this end we measured the expec- 
tation values of the Polyakov loop, the plaquette, the quark condensate ('ipifi}, and the number of 
CG matrix inversion iterations as we generated the 8 3 x 4 configurations. 

In pure gauge theory the Polyakov loop is an order parameter for the deconfinement transition: 
(P) = in the confined phase because the free energy for a single color triplet charge is infinite, 
while in the deconfined phase the test charge can be screened, so the free energy is finite and 
(P) 7^ 0. For unquenched QCD, the Polyakov loop is not an order parameter since it is nonzero 
even in the hadronic phase, but it does increase dramatically at the transition. In this work, we 
identify the thermal crossover as the place where the derivative of (Re P) is greatest. Figure ||] shows 
(Re P) versus the hopping parameter k for the seven values of fixed coupling /?; Figure [| shows 
only the runs where the crossover is at the lowest 3 values of Mps/My. Although the crossover 
becomes steeper at stronger coupling, there is no evidence of a first order transition. 
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Figure 3: Polyakov loop vs. hopping parameter - all /3's. 
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Figure 4: Polyakov loop vs. hopping parameter - for the 3 lowest /3's. 



In continuum QCD with massless quarks, one expects to see a restoration of the spontaneously 
broken chiral symmetry at high temperatures. The order parameter for this transition is the chiral 
condensate (i/>V). Since Wilson fermions break chiral symmetry explicitly, the meaning of (tpi/j) at 
k 7^ k c is not so clear. Besides the usual multiplicative renormalization one must make a subtraction 
to compensate for the additive renormalization of the quark mass. A properly subtracted {^^f) can 
be defined through an axial vector Ward identity [p8f| . However, since our study did not include 
calculation of screening propagators, we can only look at the unrenormalized ('tpip). In spite of 
these problems, Figure |5| shows a drop in {^ip) at the crossover identified by (Re P). 

Since we use the plaquette (in the space-space planes) to self-consistently tune uq, we must 
ensure that it varies smoothly across the thermal crossover. Figures || and |7| show that this is the 
case. In fact, the plaquette on the zero temperature lattices agrees within errors with the plaquette 
at finite temperatures on the confined side of the crossover. The dashed vertical lines in those 
figures simply mark the location of the crossover, kt{(3)- The large errors on the deconfined side 
are due to the smaller sample sizes where running at lower quark mass is expensive. 

As the thermal crossover line Kt{0) approaches the critical line k c ((3), the number of iterations 
needed to invert the fermion matrix per time step, Nn eT , peaks at the thermal crossover. The reason 
is that as one approaches kt(P) from the confined side (varying k with (5 fixed) the zero modes at 
k c {/3) become more influential, while there are no zero modes in the deconfined phase. Figure || 
shows the peaks in iVi ter are at the same locations as the crossovers indicated by the Polyakov loop. 
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Figure 6: Space-space plaquette vs. hopping parameter - all (3's. The shaded octagons mark the 
zero temperature values. 
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Figure 7: Space-space plaquette vs. hopping parameter - for the 3 lowest /3's. The shaded octagons 
mark the zero temperature values. The dashed lines indicate our determination of kt(P) and 
emphasize the agreement between the space-like plaquettes measured on zero temperature and 
finite temperature lattices at the crossover. 
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Figure 8: Conjugate gradient matrix inversions vs. hopping parameter. 
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4.2 Spectrum 



For a number of reasons, it is useful to evaluate some zero temperature quantities at the parameters 
of our thermodynamic simulations. The light hadron spectrum is essential in determining the chiral 
limit for Wilson- like fermions. 

The spectroscopy was an entirely straightforward lattice computation which used Gaussian- 
smeared source wavefunctions and point-like sink wavefunctions. We performed correlated fits to 
a single exponential and selected the best fits based on a combination of smallest chi-squared per 
degree of freedom and largest confidence level. Propagators are separated by 10 HMC trajectories. 

Our calculations of the hadron spectrum for our zero temperature simulations are summarized 
in Table [l], and the phase diagram (Fig. ^) illustrates the location of the T = runs with respect 
to the thermal crossover and critical lines, kt(P) and k c (P) respectively. An anomaly in Table |] is 
the small data set for (3 = 6.4, k = 0.1475. Naturally we would prefer to have more configurations 
with which to compute hadron correlators. Unfortunately the cost of running at those parameters 
is high. 



Table 1: Masses of the pseudoscalar and vector mesons and the nucleon on an 8 3 x 16 lattice. 
Points marked by as asterisk lie on the Nt = 4 crossover. The # column lists the number of 
propagators used to compute spectroscopy. 



p 


K 


Uq 


# 


aMpg 


aMy 


aM N 


Mps/My 


M N /My 


6.40 


0.145 


0.826 


81 


0.931(4) 


1.351(18) 


2.14(3) 


0.689(10) 


1.58(3) 


*6.40 


0.1475 


0.828 


30 


0.664(8) 


1.26(6) 


1.63(11) 


0.527(26) 


1.29(11) 


6.60 


0.140 


0.834 


64 


1.173(4) 


1.481(9) 


2.30(3) 


0.792(6) 


1.55(2) 


*6.60 


0.143 


0.841 


144 


0.927(4) 


1.280(8) 


1.958(18) 


0.724(6) 


1.530(15) 


6.60 


0.146 


0.855 


40 


0.468(15) 


1.04(13) 


1.34(5) 


0.45(6) 


1.29(17) 


6.80 


0.1325 


0.842 


79 


1.494(3) 


1.700(7) 


2.651(13) 


0.879(4) 


1.56(1) 


*6.80 


0.137 


0.849 


120 


1.187(3) 


1.421(6) 


2.190(10) 


0.835(4) 


1.541(10) 


6.80 


0.140 


0.857 


43 


0.885(8) 


1.182(16) 


1.75(4) 


0.749(12) 


1.48(4) 


*7.20 


0.118 


0.864 


143 


1.915(3) 


1.994(3) 


3.110(5) 


0.960(2) 


1.560(3) 


*7.30 


0.114 


0.8695 


30 


2.043(4) 


2.106(5) 


3.297(11) 


0.970(3) 


1.614(6) 



While in Figure |3] we do not see the same first-order jump in (Re P) that we did with the 
standard Wilson actions, we would like to make the comparison more convincing. After all, we 
cannot know a priori the relation between the bare parameters for the standard action (/3\\r,Kw) 
and those for the improved action (f3j, ki); it could happen that a small change in Kyy corresponds 
to a much larger change in the quark mass than does a similar change in kj, giving us the illusion 
that the crossover is broader for the improved action. 

Using measurements of the pseudoscalar mass near the crossover region, we can interpolate in 
order to estimate (aAfps) 2 as a function of 1/k for both actions. In particular, we look at the 
thermal crossover for improved and unimproved Nt = 4 Wilson actions at comparable Mps/My. 
Data from Ref. [ 30 1 suggests that we compare Pw = 4.94, kw = 0.18 where Afpg/Afy = 0.836(5) 
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Figure 9: Phase diagram for the Symanzik-improved action. Octagons represent the N t = 4 
thermal crossover, and diamonds indicate estimates of the locations of vanishing pion mass. Zero 
temperature simulations were performed at the crosses. Dashed and dotted lines are merely to 
guide the eye. 
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Figure 10: Pseudoscalar/ vector meson mass ratio vs. the lattice pseudoscalar mass squared. The 
diamonds correspond to (3\y = 4.9 simulations with the unimproved Wilson action, and the squares 
denote our (iiw = 6.8 simulations with the improved action. The solid lines indicate the region of 
the thermal crossover. 
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Figure 11: Polyakov loop as a function of (aMpg) 2 for fixed j3. Squares: (3 = 6.8 improved Wilson 
fermions, Diamonds: = 4.9 unimproved Wilson fermions. Both have similar Mps/My at the 
crossover. 



to our Pjw = 6.80, kjw = 0.137 where Mps/My = 0.835(4). The data we use for comparison to 
ours are the unimproved 0w = 4.9 data provided in Refs. fl], ^TJ. We interpolated the pseudoscalar 
meson mass using a linear least squares fit to a quadratic in 1/k around the crossover region for 
the unimproved Wilson data. Due to the smaller data sample, we interpolated (aMps) 2 linearly in 
1/k for our improved Wilson data. Figure 10 is a plot of Mps/My vs. (aMps) 2 for both actions. 
The crossover occurs at similar Mps/My as indicated by the solid lines in the graph. 

Having eliminated bare quantities, we can plot the thermodynamic observables against the 
pseudoscalar mass squared. Figure 11 demonstrates that the crossover is indeed smoother for 
improved action than the Wilson action. 



4.3 Heavy Quark Potential 

In addition to computing hadronic masses, we used Wilson loop data to measure the heavy (or 
static) quark potential V(f): 

V(r) = - lim - lnW(r,t). (8) 

t^oo t 

A standard ansatz for the form of the potential is 

V(r) = V + ar-~-f(G L (r)-± : \ (9) 
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Table 2: Fits to the heavy quark potential along the Nt 
fermions). 



4 crossover (n/ = 2 improved Wilson 



p 


K 


# 


^min '"max 


t 


aV 


a 2 a 


e 


/ 


r /a 


6.40 


0.1475 


30 


1.41-4.47 


2 


1.0(3) 


0.41(8) 


0.7(2) 


5.6(6) 


1.52(4) 


6.60 


0.1430 


108 


1.41-6.93 


2 


0.65(9) 


0.42(3) 


0.34(8) 


3.21(24) 


1.77(2) 


6.80 


0.1370 


95 


1.41-6.93 


2 


0.70(6) 


0.346(15) 


0.38(5) 


2.45(18) 


1.913(13) 


7.20 


0.1180 


117 


1.41-5.66 


2 


0.65(2) 


0.253(6) 


0.33(2) 


1.06(8) 


2.287(13) 



where Vq, a, e, and / are fit parameters, and Gl is the lattice Coulomb potential. In practice, this 
fit is performed for a fixed t; that is, the potential is estimated through an effective potential, 



Vt(f) 



In 



W(r, t + 1) 



W(r, t) 



such that 



W(r,t) ~ exp(-V t (r)t) 
The parameter ro is defined to be the length such that 

r 2 F(r ) = 1.65, 



with 



F(r) 



dV 

dr ' 



(10) 

(11) 
(12) 
(13) 



which corresponds to r$ = 0.49 fm from potential models. Sommer showed this to be a useful 
quantity with which to set the lattice scale [32|. In this work we calculated the force by taking 
numerical differences of the potential. Our analysis proceeds as in Ref. [33]. Errors are estimated 
by bootstrapping the data, and occasionally increased to account for differences in the choice of t. 
We present our fits to the potential for the zero temperature simulations along the JVj = 4 crossover 
with improved Wilson fermions in Table ^. 

For comparison, we performed the same calculation with two flavors of Kogut-Susskind fermions 
for three parameter sets along the iVj = 4 crossover. Our fits are given in Table |3| Meson masses 



were taken from Table 1 of Ref. [34]. In addition, we measured the potential at one point along the 
Nt = 6 KS crossover. That fit also appears in Table [3|. The generation and spectroscopy of those 
configurations are discussed in Ref. [|35||. 



4.4 Scaling tests 



In Sections |4.l| and [4^ we showed that thermodynamics with the improved action does not have 
the same artificial first order behavior that unimproved Wilson thermodynamics does. However, 
in order to make physical predictions which can be compared with results from Kogut-Susskind 
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Table 3: Fits to the heavy quark potential along the N 
fermions). 



4(*6) crossovers (n/ = 2 Kogut-Susskind 



p 


am q 


# 


^min ''max 


t 


aV 


a 2 a 


e 


/ 


r /a 


5.2875 


0.025 


55 


1.41-6.93 


2 


0.80(10) 


0.30(3) 


0.46(10) 


1.46(20) 


1.99(4) 


5.3200 


0.050 


67 


2.24-6.93 


2 


0.68(22) 


0.29(4) 


0.2(3) 


5.7(1.1) 


2.17(11) 


5.3750 


0.100 


90 


1.00-5.66 


3 


0.62(8) 


0.288(23) 


0.26(7) 


0.56(12) 


2.20(4) 


*5.415 


0.0125 


280 


2.24-6.71 


3 


0.76(2) 


0.130(5) 


0.36(2) 


1.0(2) 


3.14(5) 



thermodynamics, we must make use of the spectrum and potential computations described in the 
preceding two sections. 

In Figure [l^ we plot the ratio T c /My as a function of the pseudoscalar /vector meson mass 
ratio Mpg/My. Extrapolation to the physical pion/rho mass ratio is necessary in order to make 
a prediction for T c . The fact that T c /My is independent of Nt for the Kogut-Susskind action 
leads one to believe that this quantity is scaling at lattice spacing a = 1/(4T C ). Clearly, this 
statement is not true for the unimproved Wilson action. The N = 4 unimproved Wilson points 
show a large dependence on the quark mass, and disagree significantly with the corresponding 
N = 6 points at Mps/My < 0.8. In addition, since T c /My is consistently lower for the improved 
action than the unimproved action at equal lattice spacing, the discretization errors in the latter 
must be appreciable. The improved Wilson point at Mps/My = 0.53 appears to have some slight 
agreement with the Kogut-Susskind data, but with a large error. Finally, we remark that one 
expects T c /My —* in the infinite quark mass limit since the vector meson mass diverges there, 
so ultimately we want to simulate at as small Mps/My as possible in order to extrapolate to 
M^/Mp = 0.18 reliably. 

In order to look at T c scaled by quantities which are nominally independent of the quark mass, 



we use yfa and ro from our potential fits mentioned in Sec. 4.3 . In Figure [L3] the rise in T c j ' yfd and 
roT c as Mps/My — > 1 is presumably due to T c approaching the pure SU(3) transition temperature 
as the quarks decouple. The Nt = 4 quenched T c j^fa from Ref. |36[ appears as an arrow in Fig. 13 
and supports this presumption. The disagreement between the Kogut-Susskind and improved 
Wilson actions is more apparent in Fig. 13 than in Fig. 12. The error in ^fa is large, but both 
T c /\/(T and tqT c are lower for our improved action than for the KS action. In fact, the small error 
in ro reveals the presence of quark mass dependences even at Mps/My = 0.53. 



The quark mass effect can be identified further in the scaling plot of T c j ' yfa vs. a^fa (Fig. 14). 
Since a = 1/(4T C ) for all of the Nt = 4 data, the spread in ay/a for the improved action is caused by 
the increase in the deconfinement temperature as the quarks become infinitely heavy. One should 
contrast to this the observation that the 3 Nt = 4 KS points lie on top of each other. The higher 
N points for both KS and quenched actions show their relative independence on lattice spacing. 
The conclusion one should draw from Figs. 13 and [l4| is that in the case of the improved Wilson 
data, any attempt at extrapolation to physical quark mass is premature. 

In Figure 15 we plot tq^/u vs. a/rQ. While ro and y/a scale together within the error bars, the 
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Figure 12: Critical temperature divided by vector meson mass vs. pseudoscalar/ vector meson 
mass ratio. Our data for Nt = 4 improved Wilson actions are the octagons. Diamomds: iVf = 4 
unimproved Wilson; Square: = 6 unimproved Wilson; Fancy diamonds and squares: Nt = 4, 6 
Kogut-Susskind (KS), respectively M. 
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Figure 13: Critical temperature scaled by (a) the square root of the string tension and (b) the 
inverse Sommer parameter vs. pseudoscalar /vector meson mass ratio. Octagons: Nt = 4 improved 
Wilson; Diamonds: Nt = 4 KS; Square: Nt = 6 KS. The arrow in (a) shows the Nt = 4 quenched 
Tj^Ja from Ref. |§. 
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spacing in units of l/yfa. Octagons: JV^ = 4 improved Wilson; Diamonds: N% = 4 KS; Square: 



N t = 6 KS; Crosses: Quenched SU(3) from Ref. |36|]. 
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Figure 15: The dimensionless quantity r^^/a vs. the lattice spacing in units of r$. Octagons: Nf = 4 
improved Wilson, Diamonds: Nt = 4 KS. 
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Nt = 4 improved Wilson, Diamonds: Nt = 4 KS, Square: Nt = 6 KS. 



variation in a/ro for the improved action along the crossover is another manifestation of the scaling 
violations between the critical temperature and the gluonic observables. A plot against a^fa looks 
qualitatively the same, but with larger errors. 



A graph of the vector meson mass times ro (Fig. 16) shows nice behavior for the Kogut- 
Susskind simulations, disagreement between KS and clover, and the rise in My toward infinity at 
large Mps/My. Again, we do not show My/y/a vs. Mpg/My since it is qualitatively the same, 
but with larger error bars. If we were so bold as to argue that the improved Wilson data could be 
extrapolated to physical M n /M p using the points Mps/My < 0.8, then we would conclude My / y/a 
for our action is less than for the Kogut-Susskind action. This would not be too surprising given 
similar trends in scaling violations for quenched QCD spectroscopy as presented in Figure 2 of 
Ref. jilfl ) for example. At finite lattice spacing, Mp/y/a computed with KS valence quarks lie 
above the a = extrapolation, while unimproved Wilson quark calculations give a value less than 
the continuum number. The addition of the clover term significantly reduces this scaling violation; 
however, the lattice value of M p j yHj still lies below its continuum value. Of course, in the absence 
of clear scaling between My and ^Ja (and ro), such arguments in this work are speculative. 



5 Conclusions 



This is the first large scale simulation of unquenched QCD with improved Wilson fermions of 
which we are aware. We find that the Symanzik improvement program, at this level, fulfills its 
promise in that a serious lattice artifact, the spurious first-order transition at intermediate hopping 



17 



parameters, has been removed. The thermal crossover does become progressively steeper as one 
decreases the quark mass, but it is smooth in the sense that the Polyakov loop and the plaquette 
are single-valued for all (/?, k) at which we computed. 

However, improvement at this order is no panacea. It is still very costly to invert the fermion 
matrix near and below Mpg/My ~ 0.5. Since the critical temperature and the vector meson mass 
show a significant dependence on the quark mass, extrapolation to M n /M p is not trustworthy. 
Furthermore, disagreement is evident in T c /^/a and Myro between our improved Wilson action 
and the unimproved Kogut-Susskind formulation, even at comparable Mps/My. 

One cannot yet use this disagreement to cast doubt on the Kogut-Susskind results because the 
scaling behavior of the improved Wilson action has not been sufficiently tested. Simulations with 
smaller lattice spacing, perhaps Nt = 6, would give a more concrete picture of the extent of scaling 
violations in this action. 

Before beginning such an expensive undertaking, however, let us speculate as to the short- 
comings of the present action in the context of thermodynamics. In the high temperature phase, 
thermodynamic quantities are dominated by high momentum contributions. Therefore, one must 
not only improve the effects of the finite lattice spacing, but also the dispersion relation at all mo- 
menta. Although the gauge action we used has a dispersion relation closer to the continuum than 
the plaquette action, the clover term does not change the fermionic dispersion relation from that 
of the unimproved Wilson action. The work with an improved gauge action but standard Wilson 
fermions by Ref. j| shows improvement similar to ours: viz. removal of the jump-discontinuity in 
the Polyakov loop. A detailed comparison of the critical temperature from their action versus ours 
and the standard Wilson and KS actions remains to be made. 

Therefore, it is plausible that improvement of the gauge action is responsible for the removal of 
the artificial first-order behavior at intermediate values of the Wilson hopping parameter. However, 
improvement of the fermion action probably plays a role in the closer agreement to the KS results 
for T c /Mv, as was found in studies of quenched spectroscopy. Persistent quark mass dependence 
and apparent disagreement between our results and KS results for T c scaled by quark potential 
parameters indicate that further improvement in the fermionic sector is warranted. One might 
consider using Wilson-type fermions with an improved dispersion relation in the next large scale 
thermodynamics study. 
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